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摘 要 本 文 对 一 亚 音速 圆 口 喷 管内 壁面 施加 讶 孔 ， 使 其 壁面 多 孔 化 ， 采 用 大 涡 模拟 结合 FW-H 方法 研究 了 马赫 数 
M=0.6 时 喷 管 外 的 非 稳 态 流 场 与 远 场 噪声 . 研究 结果 显示 , 这 种 内 壁 小 孔 结构 对 喷 管 外 流 场 的 满 流 脉动 有 一 定 的 削弱 作 
用 . 时空 相 关系 数 的 对 比 表 明 ， 内 壁 小 孔 结构 同时 减弱 了 脉动 量 各 自在 时 间 和 空间 上 的 相关 性 ,削弱 了 大 尺度 的 涡 结构 。 


通过 对 涡 强 、 膨 胀 率 和 Q 涡 识别 准则 的 可 视 化 分 析 , 研究 了 剪 切 


层 中 涡 结构 在 喷 流 轴 向 的 产生 、 发展 与 演化 过 程 , 结果 表 


明 内 壁 小 孔 结构 抑制 了 涡 配 对 现象 , 加 强 了 流动 挨 混 ， 从 而 加 速 了 剪 切 层 的 发 展 . 


关键 词 喷 管 ; 多 孔 内 壁 ; 大 涡 模 拟 ; 非 稳 态 流 场 ; 气动 噪声 
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sonic nozzle is changed to be porous by adding 


In this paper, the inner wall 
blind holes. Large eddy simulation 5 with the FW-H method are employed to study the 


unsteady flow field and far-field noi 
the porous wall of the nozzle h 


iation at Mach number M=0.6. The results show that 
effects on reducing the turbulence intensity. Comparisons of 


the turbulence correlation rr reveal that both the temporal and spatial correlations of the 
suppressed by the porous inner wall. Large scale turbulence structures 


turbulent fluctuations arekals 
are weakened. The ge ion, development and evolution of vortices i 


by visualization and analysis of the vorticity, the dilatation and th 


the shear layer are studied 
-criterion. The results show 


that the pairing of vortices is suppressed and the flow mixing .is Gibinced by the porous wall, and as 


a result, the development of the shear layer is accelerated. N 


Key words nozzle; porous inner wall; large eddy sim 


WV 


0 前 言 


二 
大 的 影 


压 级 可 以 达到 150~160 dB 四 ,对 环境 造成 极 
响 。 在 过 去 的 半 个 世纪 里 ， 大 涵 道 比 发 动机 技术 的 
发 展 和 广泛 应 用 ， 发 动机 喷 流 噪声 已 经 得 到 了 一 定 
的 控制 。 随 着 涵 道 比 的 增 大 或 已 接近 技术 极限 趾 ， 
喷 流 噪 声 的 控制 又 重 回 了 喷 管 本 体 ， 方 法 上 大 致 包 
括 主动 方法 和 被 动 方法 . 

使 用 多 孔 结 构 对 声 源流 场 和 气动 声场 进行 控制 ， 
一 直 是 降 噪 技术 中 的 重要 方法 . 国外 将 多 和 孔 结 构 直 
接应 用 于 声 源流 场 的 报道 , 最 早 见 于 1984 年 Thiede 
等 冉 采用 多 孔 材 料 和 多 孔 结 构 改 善 咽 型 气动 性 能 
的 研究 。 类 似 地 ，Revell 等 向 试验 发 现 多 孔 处 理 
的 副 翼 侧 边 可 以 减 小 对 机 身 噪声 的 贡献 。 这 一 类 采 
收 稿 日 期 ，2016-12-27; 修订 日 期 2017-04-22 
基金 项 目 : 国家 自然 科学 基金 项 目 (No.51576067) 


作者 简介 : 林 健 (1990- ), 男 , 研究 生 , 主要 从 事 气 动 声学 研究 . 


So. unsteady field; aeroacoustic noise 


用 多 孔 材 料 结构 减弱 声 源流 场 的 应 用 还 可 见 诸 文 献 
[5]~[7] 的 报道 。 总 体 而 言 ， 虽 然 多 孔 处 理 已 被 证 明 
是 流动 与 噪声 被 动 控制 的 一 种 有 效 措施 9 , 但 是 目 
前 对 多 和 孔 流 动 控制 技术 的 研究 还 不 够 广泛 ， 对 其 机 
理 的 认识 也 不 够 清晰 . 在 喷 流 流动 以 及 噪声 控制 问 
题 中 ， 采 用 多 孔 材 料 或 多 孔 处 理 方法 的 研究 ， 无 论 
是 在 实验 测量 或 数值 模拟 分 析 方面 ， 目 前 国内 外 尚 
未 有 公开 发 表 的 文献 报道 

本 文采 用 大 涡 模 拟 与 FW-H 声 比 拟 方法 相 结 合 ， 
在 马赫 数 M = 0.6 时 , 对 三 种 不 同 内 流 道 的 喷 管 分 
别 进行 流 场 和 噪声 计算 , 对 比分 析 内 壁 小 孔 结 构 对 
喷 流 流 场 的 时 均 速度 和 脉动 量 w、v 及 p' 分 布 的 影 
响 。 在 此 基础 上 ， 进 一 步 研究 近 喷 管 出 口 剪 切 层 内 
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uw、v/'、p/' 在 时 间 和 空间 的 相关 性 , 分 析 其 分 布 与 发 
展 的 特点 。 从 锅 强 、 膨 胀 率 以 及 Q 判 据 的 可 视 化 分 
析 , 研究 该 多 孔 结 构 对 旋涡 在 喷 流 轴 向 的 发 生 、 发 
展 与 演化 过 程 的 影响 .和 远 场 噪声 的 频谱 分 析 则 通过 
声 比拟 结果 完成 ,目的 是 初步 检验 小 孔 结 构 对 辐射 
噪声 的 抑制 作用 ， 
1 数值 计算 方法 
1.1 物理 模型 及 网 格 划分 

本 文 计 算 对 象 是 三 个 有 出 口 延 伸 段 的 渐 缩 喷 管 ， 
其 中 基础 喷 管 I 的 出 口 延 伸 段 为 光滑 内 壁 , 喷 管 了 I 和 
亚 的 延伸 段 内 壁 则 有 盲 孔 ,图 1 是 喷 管 几何 模型 及 
育 孔 局 部 放大 示意 图 。 喷 管 出 口 直径 D; = 10 mm， 
长 度 工 = 170 mm。 喷 管 I 和 耳 内 壁 小 孔 参 数列 于 
表 1 


表 1 小 孔 参 数 
Table 1 Geometry parameters of the blind holes 


喷 管 孔径/ 孔 深 / 孔 心 距 / 周 向 轴 和 疝 

编号 mm mm mm 个 数 排 AN 
工 0.5 0.5 1.5 4 
III 1 1 2 4 


图 1 喷 管 几何 模型 


Fig. 1 Geometry model of nozzles 


计算 域 如 图 2 所 示 ， 其 中 喷 管 I 和 III 的 求 
解 域 还 包括 所 有 讶 和 孔 。 坐标 原点 位 于 喷 管 出 口 圆心 ， 
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YN 4 二 一 Pr 


为 280 万 、286 万 和 290 万 , 为 保证 计算 结果 的 可 对 
比 性 ， 除 内 部 网 格外 ， 本 文保 障 三 个 喷 管 外 部 的 网 
格 划 分 完全 一 致 。 
1.2 控制 方程 与 计算 方法 

为 对 喷 管 中 的 可 压缩 流动 进行 大 涡 模 拟 ， 采 用 
经 Favre 滤波 的 N-S 方程 上 0]， 


Op Opi 
一 1 
ee Op O06i; 
_ Oy 0 (5, 0;) 
Oz; Oz; < 全 
O(pEFPUWi/2) 0 ,~ ,~ 00 
= (Gijiii) 十 Oe 十 
+ Pi/2+D)i] _RHS， 


CN OX; 
语 顶 部 符号 “ ”表示 空间 滤波 , 而 “~” 表示 Favre 
加 权 滤 波 变量 ; p、w、p 和 e 分 别 表示 密度 、 速 

度 分 量 、 压 力 和 内 能 . 式 (2) 中 75j 为 亚 格子 应 力 ， 


Ti = P (Ui — ill;) (4) 


5 为 过 滤 后 的 苛 性 应 力 张 量 ; 5 是 以 过 滤 速 度 和 温度 
为 参数 的 备 性 应 力 张 量 项 ， 


8200 
旬 和 Pr 分 网 为 热 通 量 和 普 妆 特 常数 


XR 而 (了 ) Cn 而 


距 管 外 长 度 为 60Pi, 上 游 宽度 为 20D, ,下游 宽 度 2XO” RHSs 由 7 项 构成 , 具体 参见 文献 [11]. 本 文 亚 格子 


40D;; 声 源 面 宽度 为 10D; 至 30D;, 长 度 为 58 多、 


图 2 计算 域 


Fig. 2 Computational domain 


采用 结构 化 网 格 , 剪 切 层 区 域 进行 加 密 处 理 , 距 
管内 部 则 采用 O 型 网 格 . 喷 管 I~III 的 网 格 数 分 别 


模型 选用 动态 Smagorinsky-Lilly 模型 .为 使 方程 组 
封闭 补充 状态 方程 : 


j= RPT (7) 


式 中 R=287.1 J-kg-!.K-1. 

本 文 使 用 商业 软件 FLUENT 进行 模拟 , 计算 分 
两 步 进行 . 首先 对 喷 流 流 场 进 行 数 值 仿 真 ; 其 次 提 
取 声 源 面 上 的 流 场 数据 ,通过 FW-H 方程 计算 喷 流 
噪声 场 . 流 场 计算 选用 Density-Based 求解 器 ,流体 
介质 设置 为 理想 气体 ， 务 度 采 用 Sutherland 模型 , 
采用 RNG k 一 e 消 流 模型 进行 稳 态 计算 ,整体 质量 
流量 误差 小 于 1% 时 认为 稳 态 达到 收敛 。 然后 开始 瞬 
态 计算 , 为 充分 求解 涡 脱落 过 程 以 及 保证 CFL 数 小 
于 1, 设置 非 稳 态 时 间 步 长 ttep 为 5x10-"7 s, 非 稳 态 
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表 2 喷 管 模拟 结果 的 对 比 


Table 2 Comparison of simulation results of nozzles 


喷 管 编号 理想 值 
(Uj,max)#/ (m/s) 204 
(Uji)tm/ (m/s) 204 
(Pi,total)tm/ Pa 128706.9 
(pj,total) im /Pa 128706.9 
(Vptotal)tm /Pa 0 
(Vptetel) em £ (pi,total) tm 0 
(ri)tm/ (Kg/s) 0.0193 
(rj)tm/ (KE/s) 0.0193 
Th/N 3.57 


计算 时 间 达 到 tU;/D; > 200 时 认为 流动 状态 已 经 可 
以 排除 稳 态 初 场 的 影响 即 达 到 统计 学 稳定 ， 此 后 开 
始 声场 计算 并 进行 流 场 数据 采样 统计 ， 采 样 时 间 持 
续 47/Dj > 204。 
1.3 边界 条 件 

为 保证 三 种 喷 管 模型 出 口 马赫 数 都 为 0.6, 边界 


条 件 统 一 给 定 为 : 上 游 喷 管 入 口 给 定 总 压 Dtotal oN 


128706 Pa， 总 温 Tiotal = 313.875 K, 同时 给 出 
静 压 pstatic = 101325 Pa; 喷 管 外 计算 域 入 
周 环境 界面 设置 为 无 反射 的 压力 边界 
压 ptotal = 101325 Pa， 总 温 Ttotal 
考 静 压 pstatic = 101325 Pa; 下 
的 压力 出 口 边 界 条 件 , 给 定 static = 101325 Pa， 
总 温 Tiotal = 293.15 K; 喷 管 图 壁 采用 无 滑 移 、 绝 热 
壁面 边界 条 件 . 

本 文 前 期 工作 已 通过 实验 数据 的 对 比分 析 3 ， 
验证 了 计算 模型 、 方 法 以 及 边界 条 件 的 有 效 性 ， 


2 计算 结果 及 分 析 


2.1 平均 流动 特性 比较 x 


为 方便 表述 , 对 本 文 的 符号 作出 如 下 说 明和 和 
ww 分 别 表示 轴 向 速度 和 轴 向 速度 脉动 ; v 和 wv 则 表 
示 径 向 速度 和 径 向 速度 脉动 , p 和 yp/ 则 表示 压力 和 
压力 脉动 ; % 表示 质量 流量 ; Ti 表示 推力 。(.》 表示 
变量 的 平均 值 ,其 涉及 的 下 标 + 表示 时 间 平均 , 0 表 
示 周 向 平均 , m 表示 质量 加 权 平均 . 其 余下 标 1 和 j 
分 别 表示 喷 管 的 进口 和 出 口 ,c 表示 喷 管 中 心 线 . 

表 2 给 出 了 喷 管 1、 、 三 模拟 结果 以 及 理想 气 
休 完全 膨胀 值 与 之 间 的 比较 .可 以 发 现 由 于 理想 状 
况 为 无 黏 流动 ， 忽 略 了 喷 管 内 部 的 压力 损失 ， 导 致 
晓 管 出 口 处 质量 加 权 轴 向 平均 速度 、 质 量 流量 、 推 
力 都 大 于 喷 管 1 、IE 和 亚 的 情况 ， 对 比 喷 管 工 和 号 
管 了 、 亚 ,推力 分 别 为 324N.、3.23N 和 3.24 N, 即 
加 入 内 壁 小 孔 结构 后 喷 管 推力 损失 基本 可 以 忽略 不 


喷 管 I 喷 管 II 喷 管 三 
205 204 204 
197 195 196 

128706.9 128706.9 128706.9 
126917.7 126510.1 126724.9 

1789.2 2196.8 1982 

1.39% 1.71% 1.54% 

0.0181 0.0182 0.0182 

0.0181 0.0182 0.0182 
3.24 3.23 3.24 


计 , 因此 本 文 使 用 的 多 孔 参 数 较为 合理 . 
图 3(a) 为 喷 管 出 口 轴 向 速度 沿 径 向 分 布 对 比 . 


可 以 看 出 喷 管 工 出 口 高 速 流体 分 布 范围 最 大 ;小 孔 
全 移 人 RER 了 和 四 内 部 这 的 名 站 下 页 
CS 


0.4 


0.3 


0.2 


(uc)e/U; 


x/D; 


图 3 喷 管 出 口 轴 向 速度 及 中 心 线 上 轴 向 速度 分 布 
一 , 喷 管 1; - - - -, 喷 管 II; 一 .一 . 一, 喷 管 III 
Fig. 3 Distributions of axial velocity along radius and the 
center axial lines. 一 , nozzle LI; - - - -, nozzle II; 一 :一 :一 


nozzle III 
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速 流 体 分 布 范围 都 略 有 减 小 .同时 如 图 3(b) 所 示 ， 
喷 管 I 和 左 沿 径 向 速度 梯度 减 小 ， 这 减弱 了 喷 流 向 
下 游 发 展 过 程 中 卷 入 周围 流体 介质 的 能 力 ， 使 得 主 
流 与 周围 环境 发 生 质 量 、 动 量 和 能 量 的 交换 变 得 没 
有 了 喷 管 I 那样 剧烈 ,因此 出 现 了 下 图 中 噶 管 I 、U 
的 势 流 核 长 度 不 减 反 增 的 现象 ， 
2.2 脉动 特性 比较 

图 4 给 出 了 不 同 轴 向 位 置 潜流 脉动 沿 径 向 的 分 
布 . 可 以 看 出 在 z = Di 轴 向 位 置 处 哎 管 I 和 喷 管 王 
的 轴 向 速度 脉动 w 略 大 于 喷 管 I 往 下 游 发 展 过 程 
中 两 喷 管 的 wu/ 都 得 到 了 削弱 , 到 x = 5D; 位 置 时 都 
已 减 小 低 于 喷 管 I . 径 向 速度 脉动 v 的 减弱 趋势 比 
的 减弱 趋势 更 明显 , 在 x = Di 轴 向 位 置 时 ， 喷 管 
I 和 本 的 v 就 已 降 至 低 于 喷 管 I 的 值 . 脉动 速度 分 
量 乘 积 wv 以 及 压力 脉动 p' 同样 体现 出 了 一 定 的 减 
弱 趋 势 。 
2.3 脉动 的 时 空 相关 性 分 析 

流 场 脉动 的 时 空 全 相关 系 儿 可 用 来 考 守 丰 同时 和 NA 
或 不 同位 置 两 点 脉动 基 间 的 相关 性 . 如 图 5 所 示 
于 剪 切 层 中 脉动 值 ,点 z 与 其 下 游 点 z 十 sk 
时 空 相关 系数 Ri 定义 如 下 : 


Ri (RE TT) = 


式 中 , z 为 测 点 在 流 场 中 的 轴 向 坐标 ,，: 为 两 个 测 点 
间 的 距离 ,7 为 延迟 时 间 ， 

图 6(a)~(c) 分 别 表示 喷 管 层 线 上 x = 2.5D; 处 
轴 向 脉动 速度 w、 径 向 脉动 速度 v De . 
在 时 间 序列 上 的 自 相 关 性 Ruww、Rws 和 已 mw。 \\ 
小 孔 结构 对 此 处 Ru 的 影响 其 微 ,而 对 Ro。 和 忌 RN 
的 影响 较 明 显 , 它们 各 自 的 时 间 相关 性 被 削减 . SS 


Et 


林 健 等 : 喷 管内 壁 多 孔 处 理 的 圆 口 喷 流 大 涡 模 拟 


A (c) 70(u 


4 5 
v /UF 


| 
| 
| 
3 


ee / Go 


图 4 不 同 轴 向 位 置 油 流 脉动 沿 径 向 的 分 布 . 每 幅 图 中 从 左 到 右 
的 位 置 /Di 分 别 为 1、2.5 和 5, 线 型 含义 同 图 3 


Fig. 4 Radial profiles of turbulent fluctuations at axial 


positions xz/D; =1, 2.5 and 5. Legend as in Fig. 3 


地 ,本 文 小 孔 结 构 对 压力 脉动 p 分 别 与 w、v 在 时 
间 序 列 上 的 关联 Ru 和 Rs 也 有 明显 的 削弱 作用 . 
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时 间 相关 性 的 减弱 一 定 程度 上 反映 了 流 场 混乱 程度 ”减弱 了 脉动 量 w、w 和 py/ 各 自在 空间 上 的 相关 性 ， 
的 增强 ， 说 明 内 壁 小 孔 结构 加 速 了 大 尺度 旋涡 向 小 ”再 次 从 空间 角度 印证 了 小 孔 结构 减 小 旋涡 结构 的 
涡 的 转变 . 作用 . 

图 7 表示 喷 管 展 线 上 x = 10D; 处 三 个 喷 管 的 
Rw、Ros 和 Rss 所 对 应 的 云图 . 图 中 虚线 斜率 代表 一 一 一 一 一 :下 
了 该 位 置 脉动 量 向 下 游 的 传递 速率 . 可 以 发 现 喷 管 
I 和 三 中 虚线 的 斜率 大 于 喷 管 了 中 的 对 应 值 , 可 见 。 二 一 


内 壁 小 孔 结 构 加 速 了 w、v 和 yp 向 下 游 传 递 的 速率 图 5 前 切 层 中 脉动 相关 系数 定义 的 对 应 位 置 zx 和 测 点 间 
. i : 距离 示意 图 
此 外 ， 喷 管 的 时 空 相关 系数 正 值 区 域 比 喷 管 I 和 Fig. 5 Two-point spatial-temporal correlations in the shear 


I 较 宽 了 说 明 在 同 样 的 时 间 间 隔 下 3 内 壁 小 孔 结 构 layer, obtained for locationz and using a spatial separation é€ 


3 4 5 6 7 


图 6 喷 管 唇 线 上 z = 2.5D; 处 脉动 相关 系 3 二 班 信义 同 图 ， 3 


Fig. 6 Two-point spatial-temporal correlations of fluctuations i i hear layer at x = 2.5D;. Legend as in Fig. 3 
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图 7 喷 管 层 线 上 x = 10D; 处 脉动 相关 系数 云图 
Fig. 7 Contours of Ruu, Rvv and Rpp at x = 10D; in the shear layer 
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2.4 前 切 层 厚度 分 析 

前 切 层 厚度 6 = (ro0.10 一 70.95), 其 中 rol 和 ro.95 
分 别 表 示 轴 疝 流 速 为 10%(wc); 和 95%(wc); 处 对 应 的 
半径 值 . 图 8 为 剪 切 层 厚 度 对 比 图 , 3 条 分 布 曲线 基 
本 都 呈现 出 线性 增长 的 趋势 , 这 与 剪 切 层 理 论 一 致 . 
对 比 其 斜率 ， 可 以 发 现 内 壁 小 孔 结 构 降 低 了 剪 切 层 


厚度 增长 的 速度 . 


Ou 


Ov 


Ow 


O00 bag fangu wov) 
下 和 队 Oz Ovy Oy Oz Oy Oz 
CN FLT 

人 oz Oz Ox Oz 


图 9 给 出 了 7 一 xz 平面 上 由 涡 强 |2| 和 膨胀 
率 表 征 的 瞬 态 涡 结构 图 .图 9(a)~(c) 可 以 看 出 距 
离 喷 管 出 口 约 1.5 倍 半径 范围 内 喷 管 都 存在 一 个 准 
层 流 环形 混合 层 . 这 一 准 层 流 环形 混合 层 的 发 展 主 
要 受 不 稳定 性 波 引起 的 扰动 增长 的 影响 ， 导 致 剪 切 
层 变形 以 及 旋涡 卷 起 形成 大 尺度 涡 环 结构 ， 随 后 涡 
环 向 下 游 对 流 时 产生 涡 配 对 现象 ， 它 也 是 喷 流 噪声 
产生 的 原因 之 一 ;最 终 涡 环 失去 稳定 性 而 变形 . 图 
9(a)、(b) 中 虑 组 窒 标 出 的 地 方 可 以 清晰 的 看 到 涡 
配对 这 种 非 钱 性 相互 作用 现象 , 而 在 图 9(c) 中 这 种 
现象 六 和 P 明 显 。 比较 图 9(a)~(c) 的 膨胀 率 云图 . 可 

， 喷 管 I 和 工分 布 规律 类 似 ， 在 喷 管 出 口 附 


以 
0 1 2 3 4 5 6 7 8 9 10 NE 让 传播 方法 主要 是 上 游 和 径 向 方向 ， 而 喷 管 开 


图 8 前 切 层 厚度 对 比 图 线 型 信义 同 图 3 AGS 
Fig. 8 Comparison of shear layer thickness. 
Legend as in Fig. 3 


2.5 内 壁 小 孔 结构 对 流 场 影响 的 分 析 

Powell 43 和 Howe 64 发 展 的 涡 声 理论 表明 , 声 
波 的 产生 同 流 体 中 的 旋涡 和 势 流 、 以 及 旋涡 之 间 的 
相互 作用 密切 相关 . 本 文 从 涡 强 | 如、 膨胀 率 D 和 
Q 值 这 三 个 物理 量 出 发 , 研究 多 孔 壁 面 对 流 场 中 旋 
涡 结构 和 强度 等 方面 发 挥 的 作用 ,重点 考察 前 切 层 
中 大 尺度 涡 结构 的 发 生 、 发 展 与 演化 过 程 ， 研 究 考 
察 K-H 不 稳定 性 的 发 展 与 旋涡 配对 . |Q2|, D 和 Q 的 
定义 计算 式 为 


I2| = |v x ?| (9) 


) 人 则 出 现 了 一 部 分 声波 往 下 游 传播 ， 说 明 喷 管 亚 这 组 


多 孔 参 数 已 经 对 旋涡 结构 的 改变 产生 了 较为 明显 的 
作用 . 

为 观察 消 流 涡 结构 沿 流向 的 三 维 演化 特征 以 及 
开 孔 对 涡 结 构 产生 的 影响 , 图 10 给 出 了 涡 量 在 不 同 
的 流向 位 置 > 横 截面 上 的 分 布 . 在 z/D; = 0.5 处 
可 以 看 出 轴 对 称 涡 环形 成 ， 这 与 线性 稳定 性 理论 的 
结果 一 致 4519. 进一步 的 发 展 , 旋涡 出 现 周 向 的 弯 
曲 而 开始 失去 稳定 性 ,随后 涡 环 变形 ,失去 对 称 性 ， 
形成 涡 因 结构 。 喷 管 了 和 喷 管 全 由 于 径 向 扰动 的 加 
入 ， 加 速 了 剪 切 层 的 发 展 .因此 轴 对 称 涡 环 的 特征 
不 再 明显 , 在 0.5D; 处 (图 10(d)、(g)) 涡 环 已 开始 
变形 , 而 喷 管 1 在 z/D; = 2.5 处 才 出 现 的 涡 办 结 检 
(图 10(c)) 在 喷 管 了 和 匡 中 提前 至 z/D; = 1 处 即 已 
出 现 . 
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(a) 喷 管 


ChinaXiv 合 作 期 刊 
38 着 


Xr; 
(c) 喷 管 亚 


(b) 喷 管 工 


图 9 7 一 zx 平面 上 由 涡 强 |12| 和 膨胀 率 表征 的 瞬 态 涡 结 构图 涡 强 和 膨胀 率 灰 度 云图 由 亮 到 暗 表 征 的 范围 分 别 为 (0.001,10) 和 
(一 0.001;0.001) 
Fig. 9 The instantaneous vortices shown by |2| and dilatation field in 7 一 xz plane. The gray contour levels from light to dark are 
scaled by (0.001,10) for vorticity and (一 0.001,0.001) for dilatation 


(g) 


(h) 0) 


图 10 这 人 而 和 


Fig. 10 The structtre of vortices at different x positions 


图 11 分 别 展 示 了 三 个 喷 管 模型 的 @ 准则 等 值 
面 分 布 , 其 中 @ = 2000000, 灰 度 代表 了 轴 向 流速 的 
大 小 . 从 图 11 中 可 以 看 出 在 近 喷 管 出 口 区 域 ， 剪 切 
层 扰动 涡 呈 现 周 向 的 管状 结构 ， 这 种 前 切 层 不 稳定 
性 属于 开尔文 - 效 姆 稚 兹 不 稳定 (KHI)。 随 后 前 切 
层 中 的 扰动 涡 出现 明 显 的 周 向 破碎 ,并 向 下 游 发 展 
演化 .下游 涡 结构 表现 出 三 维特 性 ， 比 如 项 链 涡 , 展 
展 向 涡 以 及 沿 流 向 细 长 的 流 线 涡 . 对比 图 11 中 局 部 
放大 图 可 以 发 现 ， 喷 管 工 和 喷 管 王 周 向 管状 结构 斥 
度 明显 小 于 喷 管 I 中 出 现 的 尺度 ,再 次 证 明了 内 壁 


小 孔 结 构 加 强 了 流动 摊 混 ， 加速 了 剪 切 层 的 发 展 
2.6 远 场 声 压 

本 文通 过 FW-H 方程 对 远 场 监测 点 的 声 压 进行 
计算 ， 得 出 不 同 监测 点 处 的 总 声 压 级 对 比 。 噪 声 监 
测 点 采用 同心 圆 面 布 置 ， 圆 心 为 喷 管 出 口中 心 ， 测 
点 距离 为 100D;, 角度 从 10° 到 150?. 图 12 是 总 声 
压 级 的 对 比 。 可 以 发 现 内 壁 小 孔 结 构 对 降 噪 有 一 定 
的 效果 ， 且 喷 管 开 的 降 噪 效 果 好 于 喷 管 I ,例如 在 
30° 处 喷 管 了 I 和 开 比 喷 管 I 的 噪声 分 别 降 低 1dB 和 
1.7 dB., 
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Streamwise Vortex 


(a) 喷 管 
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(c) 喷 管 亚 


| 据 表 征 的 瞬时 三 维 涡 结 构 及 局 部 放大 图 


Fig. 11 Instantaneous to ensional vortices represented by Q-Criterion and their local zooms 
» 


5 期 
Spanwise Vortex 
KHI 
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SPL/dB & 
图 12 监测 点 方位 示意 以 及 总 声 压 级 对 比 图 
图 中 线 型 含义 同 图 3 
Fig. 12 Positions of monitoring point and comparison of total 


sound pressure level. Legend as in Fig. 3 


为 进一步 研究 内 壁 小 孔 结构 对 声场 频谱 特性 的 
影响 , 本 文选 取 30" 方位 监测 点 进行 频谱 分 析 . 图 13 
是 该 点 1/3 倍 频谱 图 .可 以 看 出 喷 管 I 噪声 主 频 在 
2500 Hz 的 低频 部 分 . 加 入 内 壁 小 孔 结构 后 , 喷 管 I 
主 频 增 大 至 4000 Hz。 此 外 喷 管 工 和 喷 管 亚 在 1200 
Hz 以 下 的 低频 噪声 得 到 了 削弱 ， 而 1200 至 4500 
Hz 之 间 噪 声 略 有 增强 。 该 变化 的 原因 可 能 是 径 向 
扰动 加 速 大 尺度 涡 破碎 ， 相 应 地 增加 了 小 尺度 涡 的 


SPL/dB 


103 104 
Frequency/Hz 


Frequency/Hz 
图 13 30° 方位 监测 点 声 压 频 谱 图 


Fig. 13 Sound pressure spectrum at the 30° azimuthal 


monitoring point 
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3 2 论 Around a Circular Cylinder Having Porous Surface [J]|. 


本 文采 用 大 涡 模 拟 结合 FW-H 的 方法 对 具有 内 
壁 盲 孔 结构 的 喷 管 喷 流 进行 了 研究 , 得 到 如 下 结论 : 

1) 加 入 内 壁 小 孔 结构 ,由 于 沿 程 损失 和 局 部 损 
失 的 增加 ， 减 小 了 喷 管 出 口 截 面 上 高 速 流体 沿 径 向 
的 分 布 范 围 ， 同 时 沿 径 向 速度 梯度 的 减 小 ， 势 流 核 
长 度 增 加 ， 

2) 内 壁 小 孔 结构 中 流动 和 压力 的 重新 分 布 削 弱 
了 喷 管 外 流动 的 清流 脉动 

3) 喷 管 近 出 口 区 域 时 间 序 列 上 的 相关 性 
Ry，Rpp，Rup 和 Rop 因 小 孔 结构 的 影响 而 减弱 ， 
而 wv 和 p' 往 下 游 的 传递 速度 则 有 所 增 大 ， 相 应 
地 也 减弱 了 它们 各 自 的 空间 相关 人性， 

4) 喷 流 脉动 特性 的 改变 ,宏观 上 体现 为 剪 切 层 
中 涡 配 对 现象 被 抑制 和 流动 摊 混 的 加 强 ， 从 而 加 速 
了 剪 切 层 的 发 展 , 远 场 噪声 也 在 一 定 程度 上 被 抑制 。 
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